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Broadband  Plasma  Impedance  Measurements  and  Determination 

of  Plasma  Parameters 

David  D.  Blackwell**,  David  N.  Walker,  and  William  E.  Amatucci 

U.S.  Naval  Research  Laboratory ,  Plasma  Physics  Division  Washingtonf  DC 

A  small  spherical  probe  is  used  in  conjunction  with  a  network  analyzer  to  de¬ 
termine  the  impedance  of  the  probe-plasma  system  over  a  wide  frequency  range. 
Impedance  curves  are  in  good  agreement  with  accepted  circuit  models  with  plasma- 
sheath  and  electron  plasma  frequency  resonances  easily  identifiable.  Phase  measure¬ 
ments  show  clear  transitions  between  capacitive  and  inductive  modes  as  predicted  by 
the  model.  Sheath  thickness  and  absolute  electron  density  are  determined  from  the 
location  of  these  transitions.  In  addition,  much  larger  power  absorption  is  observed 
at  the  sheath  plasma  resonance  than  is  predicted  by  collional  absorption. 


I.  INTRODUCTION 


The  impedance  of  a  probe  in  a  plasma  as  a  means  of  obtaining  information  about  the 
plasma’s  characteristics  is  an  attractive  diagnostic  for  use  in  situations  where  more  con¬ 
ventional  Langmuir  probes  and  charge  collectors  are  inadequate.  In  early  experiments  by 
Harp  and  Crawford  [1]  investigating  this  technique,  the  Q  of  a  spherical  probe  immersed 
in  a  plasma  was  used  to  estimate  the  plasma-sheath  resonance  frequency.  Harp  explained 
that  a  caveat  of  this  measurement  technique  was  that  the  impedance  of  the  cable  connected 
to  the  probe  has  to  be  small  compared  to  the  impedance  of  the  sphere.  In  addition,  there 
was  no  experimental  data  given  regarding  the  plasma  frequency  measurement,  the  reso¬ 
nance  at  which  was  described  as  being  difficult  to  define.  Subsequent  authors  [2],  [3]  have 
performed  experiments  where  the  sheath  plasma  resonance  was  observed  by  changing  the 
plasma  density  while  keeping  the  frequency  of  the  signal  applied  to  the  probe  constant.  Dote 

[4]  measured  the  resonant  frequencies  of  planar  and  spherical  probes  in  a  magnetized  plasma 
by  measuring  the  DC  component  of  the  probe  current  while  sweeping  the  frequency.  More 
recently,  plasma-sheath  resonances  have  been  observed  in  a  capacitively  coupled  discharge 

[5] . 


Quantitative  measurements  of  the  impedance  itself  are  more  rare  to  find  in  the  literature 
because  of  the  difficulties/ associated  with  eliminating  stray  capacitances  and  impedances 
in  the  probe  circuit.  The  best  examples  can  be  found  in  sounding  rocket  experiments  on 
ionospheric  plasmas.  Oya  [6]  included  the  capacitance  of  the  metallic  boom  connected  to 
a  spherical  probe,  compensating  for  stray  capacitance  using  a  bridge  circuit.  Over  the  fre¬ 
quency  range  of  the  experiment,  0.5-16  MHz,  the  boom  could  be  considered  electrically  short 
enough  so  that  its  inductance  could  be  neglected.  Jensen  [7]  measured  changes  in  electron 
density  using  a  phase  locked  loop  circuit  to  track  the  frequencies  at  which  a  small  antenna 
immersed  in  a  plasma  had  an  impedance  that  was  large  and  nonreactive.  More  recently 
Steigies  et.  al.  [8]  used  a  capacitive  bridge  type  circuit  with  a  strip  antenna.  However,  all 
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of  these  experiments  were  performed  in  low  density  plasmas  such  that  u;pe/27r  is  less  than 
20  MHz.  They  also  required  specially  designed  electronics  to  overcome  circuit  effects.  In 
this  paper  we  extend  the  utility  of  this  diagnostic  method  to  much  higher  frequencies  (1 
GHz)  corresponding  to  higher  density  plasmas.  Impedance  measurements  are  made  with 
an  off-the-shelf  network  analyzer;  no  custom  parts  or  electronics  are  needed  except  for  the 
spherical  probe  itself.  Following  the  calibration  procedure  outlined  in  the  owners  manual  of 
the  network  analyzer  is  all  that  is  required  to  make  use  of  this  diagnostic. 


II.  THEORY  AND  EXPERIMENTAL  TECHNIQUE 


A.  Theoretical  Impedance 


A  spherical  capacitor  of  radius  R  surrounded  by  a  uniform,  infinite  plasma,  biased  at  the 
plasma  potential  Vp ,  will  have  a  capacitance  given  by 


C  —  AireotTR 


(1) 


with 


Cr  =  1  - 


w: 


(2) 


U}  (u>  —  ju) 

where  v  is  the  toted  electron  collision  frequency.  If  the  sphere  is  at  a  potential  different  from 
Vp,  a  sheath  will  form  between  its  surface  and  the  plasma.  The  overwhelming  majority  of 
the  sheath  impedance  is  capacitive,  so  the  total  impedance  of  the  probe  and  plasma  can  be 
written  as  if  it  were  a  capacitance  in  series  with  a  tank  circuit,  as  shown  in  figure  1.  The 
circuit  elements  are  given  by  a  sheath  capacitance  C,h,  the  vacuum  capacitance  Co  =  47t e0R, 
an  inductor  Lp  =  ui~ 2Cq1,  and  a  resistance  R  =  uLp,  with  the  impedance  given  by 


7  = 


juCth 


+ 


1 


jwCo  + 


1 


(3) 


Rp  +  juiLp 


The  sheath  capacitance  represents  the  capacitance  of  the  region  between  the  probe  surface 
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sh 
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FIG.  1:  Circuit  model  of  the  probe  plasma  impedance. 

and  the  plasma,  while  the  tank  circuit  represents  the  impedance  between  the  sheath  edge 
at  the  probe  to  the  sheath  edge  at  the  chamber  wall.  Implied  in  this  model  is  that  the 
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plasma  is  uniform  over  a  scale  length  much  larger  than  the  probe  radius.  This  circuit  model 
is  identical  to  that  outlined  in  Lieberman  [11]  for  a  parallel  plate  capacitive  discharge;  the 
only  modifications  are  that  a  spherical  geometry  is  used  and  the  sheath  capacitance  Cx , 
representing  the  sheath  between  the  vacuum  vessel  and  the  plasma,  is  neglected,  owing  to 
the  much  larger  area  of  the  space  chamber  with  respect  to  the  probe  surface  area. 

Equation  3  can  be  rearranged  using  the  definitions  of  Rp  and  Lp  and  defining  7  =  u>/u ip 
and  S  =  v /u]p  to  write 


Re  (Z) 


_L _ 8 

upCq  (7 2  —  l)2  4-  7 252 


Im(Z )  = 


Co  [(72  -  I)'  +  72*2]  +  C.h  [(72  -  1)  72  +  72<H 
u}pCoCah.  [(72  —  l)2  +  72^2 


(4) 

(5) 


The  magnitude  and  phase  derived  from  4  and  5  are  plotted  in  figure  2  for  5  =  0.01.  We  can 
see  that  there  are  two  frequencies  of  interest,  and  72,  where  the  phase  of  the  impedance 
changes  from  — 7r/2  to  7r/2,  corresponding  to  a  null  and  a  peak  in  the  magnitude.  For  small 
5 ,  we  can  solve  for  7^2  to  find  the  corresponding  points  in  frequency  space  as 


Co 


Wl 


Cah  +  Co 


W2  =  U)p 


(6) 

(7) 


If  we  assume  that  the  sheath  is  a  vacuum  region  between  the  surface  of  the  probe  and  the 
plasma  with  thickness  s  then  uji 


B.  Measurement  of  Impedance 

The  complex  impedance  Z  is  obtained  by  using  a  HP8735D  network  analyzer  to  measure 
the  reflection  coefficient  T,  given  by 

=  2j-_Zo 

with  Zo  the  internal  impedance  of  the  analyzer,  which  in  this  experiment  is  500.  Typically 
the  sphere  will  be  connected  to  a  probe  mount  with  500  coax  connected  between  it  and 
the  network  analyzer.  The  network  analyzer  will  see  a  terminal  impedance  of  this  entire 
assembly,  not  just  the  sphere  at  the  end.  The  extra  impedance  from  the  cable  and  probe 
shaft  can  be  compensated  for  calibrating  the  analyzer  to  ignore  everything  up  to  the  very 
tip  where  the  sphere  is  attached.  The  probe  should  therefore  be  constructed  such  that  the 
sphere  is  removable  and  a  short  or  50  0  load  can  be  put  in  its  place  with  no  additional 
cabling.  Figure  3  shows  how  the  probe  tip  is  constructed  from  1/4  ”  semi-rigid  coax  with 
the  outer  jacket  and  dielectric  removed  to  expose  a  short  length  at  the  tip.  One  at  a  time, 
the  operator  indicates  that  there  is  an  open  circuit  (V  =  1),  a  short  circuit  (T  =  —1),  or  a  50 
load.  The  internal  computer  then  adjusts  internal  phase  and  amplitude  filters  until  these 
values  of  T  are  visible  at  the  terminal.  Figure  4  shows  the  measured  phase  and  amplitude 
obtained  as  well  as  the  expected  theoretical  values  for  the  19  mm  diameter  sphere  used. 


(8) 
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FIG.  2:  Magnitude  (a)  and  phase  (b)  of  the  impedance  of  a  19  mm  diameter  metal  sphere  immersed 
in  a  cold  (Te  =  0.5eV)  unmagnetized  plasma  with  v/wp  =  0.01  and  a  sheath  width  of  5Aj p. 


dielectric 


FIG.  3:  Illustration  of  the  spherical  RF  probe.  The  calibration  is  done  by  removing  the  sphere 
and  attaching  a  50  Cl  load  or  short  to  the  exposed  tip  with  the  return  lead  as  short  as  possible 
and  attached  to  the  outer  jacket  of  the  coax.  An  open  circuit  is  calibrated  with  nothing  being 
attached;  this  eliminates  almost  all  of  the  contribution  from  the  capacitance  between  the  sphere 
and  the  jacket. 


C.  Experimental  Setup 

The  Space  Physics  Simulation  Chamber  at  NRL  has  been  described  in  several  earlier 
experiments  [9].  A  very  large  vacuum  vessel  (~  2-m  diameter  by  5-m  length)  is  filled  to  a 
pressure  of  p  «  10~4  Torr  Argon  and  surrounded  by  a  weak  magnetic  field  (B0  ~  3  gauss) 
such  that  <<  u,u}pe.  A  uniform,  weakly  ionized  ( n  m  107  —  109)  plasma  is  created  with  a 
large  array  of  glowing  tungsten  filaments  biased  at  -100  V.  The  electron  and  ion  temperatures 
are  typically  Te  ~  0.5  eV  and  T,  «  0.05  eV.  The  electron  density  is  measured  with  a  heated 
Langmuir  probe  mounted  in  proximity  to  the  sphere.  A  drawing  of  the  experiment  is  shown 
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FIG.  4:  Measured  magnitude  (a)  and  phase  (b)  of  the  impedance  of  the  RF  probe  in  vacuum 
compared  with  the  calculated  values  for  a  19  mm  metal  sphere. 


[SI  IS! 


ISI 
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in  figure  5. 


III.  RESULTS 


Figure  6  shows  the  amplitude  and  phase  of  the  probe  impedance  in  a  typical  plasma  with 
ne  =  108cm“3  and  a  pressure  of  3.4  X  10“4  Torr.  For  comparison  the  theoretical  calculations 
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FIG.  6:  Measured  magnitude  (a)  and  phase  (b)  of  the  impedance  of  RF  probe  immersed  in  a 
plasma  with  ne  =  2.5  X  108  cm~3  and  Tc  =  0.61  eV.  The  dashed  lines  are  computed  impedance 
curves  from  equation  3. 

are  also  shown.  There  is  a  slight  shift  in  frequency  most  likely  due  to  error  in  estimating 
the  collecting  area.  The  resonant  frequency  U72  is  plotted  versus  electron  density  inferred 
from  the  Langmuir  probe  in  figure  7.  Also  shown  is  the  plasma  frequency  derived  from  the 
formula  cupe/27r  =  8.98  x  10 3v/n7-  The  plasma  frequency  calculated  from  Langmuir  probe 
is  an  average  of  15%  lower  than  that  of  the  spherical  RF  probe  but  does  follow  the  correct 
trend. 

Figure  8  shows  the  null  frequency  oq  versus  the  resonance  frequency  u>2  as  the  electron 
density  is  increased.  For  comparison  we  also  show  u)i  =  — ,  —  •  The  sheath  thickness  s 

yrrsr 

is  the  assumed  to  be  some  number  k  Debye  lengths  \D.  The  average  value  of  s  calculated 
from  uq/u;2  gives  a  value  of  k  =  4.5,  which  is  on  the  order  of  that  estimated  by  Harp  and 
Crawford  [1].  When  a  much  larger  sphere  (R  =  7.6  cm)  was  substituted  for  the  smaller 
sphere,  the  null  frequency  uq  became  so  small  that  the  resolution  A u>  was  a  significant 
fraction  of  the  measurement.  This  is  most  likely  the  reason  for  larger  errors  seen  in  the 
plot  (figure  9(b))  of  uq  versus  o?2  corresponding  to  figure  8.  Figure  9(a)  shows  the  resonance 
frequency  u>2  is  still  a  good  approximation  of  the  plasma  frequency.  An  interesting  feature  of 
the  impedance  probe  is  that  the  reflection  coefficient  T  also  has  a  dip  at  co  =  uq ,  indicating 
a  large  deposition  of  RF  power  (figure  10).  This  dip  is  what  is  commonly  referred  to  as 
the  sheath-plasma  resonance.  The  theoretical  curve  corresponding  to  the  plasma  conditions 
is  shown  on  the  second  y-axis.  Although  the  location  of  the  resonance  is  the  same,  the 
experimental  dip  is  much  deeper  and  broader  than  what  theory  predicts,  requiring  at  least 
an  order  of  magnitude  increase  in  our  estimate  of  the  collision  frequency. 
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FIG.  7:  Electron  density  vs.  impedance  resonance  frequency.  The  dashed  line  is  the  theoretical 
value  calculated  from  the  formula  wpe/2n  =  8.98  X  10 3y/n7 


FIG.  8:  The  plasma  sheath  resonance  frequency  (uq)  vs.  the  impedance  resonance  frequency  (w2)- 
The  dashed  line  is  the  theoretical  value  inferred  from  probe  measurements  and  assuming  a  sheath 
width  of  s  =  4.5Af) 
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FIG.  9:  Relationships  between  /i,  /2,  and  ne  for  a  large  (76  mm)  diameter  sphere. 


FIG.  10:  Magnitude  of  the  reflection  coefficient  vs  frequency  for  the  plasma  conditions  of  figure  6. 
Top:  Theoretical  curve;  bottom:  measured  curve. 


IV.  CONCLUSIONS 

We  have  demonstrated  the  use  of  the  plasma  impedance  probe  at  high  frequencies  using 
a  commercial  network  analyzer.  The  probe  is  simple  to  construct  and  can  be  used  remotely 
with  a  long  transmission  line  length  between  the  analyzer  and  the  probe.  It  is  ideal  as  a 
diagnostic  in  a  plasma  with  large  amounts  of  contaminants  or  reactive  chemicals  that  would 
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interfere  with  the  collecting  surface  of  a  Langmuir  probe.  In  this  regard  it  is  similar  to 
the  plasma  absorption  probe  built  by  Kokura  [10];  however,  the  theory  and  data  analysis 
is  less  complicated.  The  density  indicated  by  the  resonance  in  the  impedance,  given  when 
iv  =  u)pe,  is  on  average  within  25%  of  the  density  measured  using  a  heated  Langmuir  probe. 
The  location  of  the  impedance  null,  refered  to  commonly  as  the  plasma  sheath  resonance, 
agrees  with  the  model  proposed  by  Harp  and  Crawford  [1].  This  gives  a  sheath  width  of 
approximately  4.5  Debye  lengths,  which  is  also  in  agreement  with  earlier  results.  The  large 
amount  of  energy  deposited  at  this  resonance  indicates  a  collsion  frequency  approximately 
an  order  of  magnitude  higher  than  theoretical  estimates. 
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